Hawking radiation from laser filaments? 
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Belgiorno et al have reported on experiments aiming at the detection of (the analogue of) Hawking 
radiation using laser filaments [F. Belgiorno et al, Phys. Rev. Lett. 105, 203901 (2010)]. They sent 
intense focused Bessel pulses into a non-linear dielectric medium in order to change its refractive 
index via the Kerr effect. We study this scenario with the parameters of their experiment and 
find that the effective metric does not quite correspond to black hole evaporation but obeys more 
similarities to cosmological particle creation - even though this effect can probably not quite explain 
the observations. Modelling Hawking radiation would require much stronger pulses. 



PACS numbers: 04.62.+V, 42.65.Hw, 42.65.Ro, 04.70.Dy. 



I. INTRODUCTION 

The idea of finding analogies to gravity in labora- 
tory systems is quite old. For example, already in 1923, 
W. Gordon used an effective metric in order to describe 
light propagation in media Later it has been real- 
ized [2[ that these analogies could be employed to model 
quantum particle creation phenomena such as Hawking 
radiation 0, 01 • The Hawking effect is the thermal emis- 
sion of radiation from horizons - cither of black holes in 
general relativity, or in analogue systems, where such a 
horizon is that of the effective space-time metric defined 
by the wave equation for light (or other waves) propagat- 
ing in the medium. An analogue of a black hole horizon 
can occur where the velocity of the medium exceeds the 
speed of light (or sound) within the medium, thereby 
dragging along the photons (or phonons). 

By now there are many proposals for such analogue 
systems based on various laboratory systems, for instance 
sound propagation in flowing fluids, see, e.g., A 
black hole analogue for photons has been proposed in Q 
employing the phenomenon of slow light - for a discussion 
of the prospects and difficulties of slow light, see, e.g., 
[l0l - [l3| . Partly motivated by these problems, black hole 
analogues based on the propagation of light in ordinary 
dielectrics were discussed in [ijj [3 ■ An important idea 
was the realization that the motion of the medium can 
be replaced by the motion of a pulse thr ough the medium 
which changes the local light speed, cf. [l^ ■ 

Recently, Belgiorno et al have argued that they 
have seen Hawkirig radiation from such an analog sys- 
tem, see also [3, [l^ ■ Let us review their claim. They 
direct an intense pulse of radiation into a piece of sil- 
ica glass, such that in a very thin region of finite size 
the Kerr effect increases the refractive index of the glass 
n = uq + Sn. At a frequency lu separated from the frc- 
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quency of that incoming radiation, the phase velocity 
fph(w) of light in the medium is such that outside the 
pulse Wph(w,?T.o) is (slightly) greater than the speed of 
the pulse - while within that region, the phase velocity 
Vph{u!,no + Sn) is slightly less than the pulse speed (due 
to the increased refractive index no+Sn). The point were 
the phase velocity (at that particular frequency) equals 
the velocity of the pulse is claimed to be a "phase hori- 
zon" . Looking at a direction perpendicular to the pulse, 
they see a slight excess of photons (about one photon ev- 
ery ten pulses) at this frequency where their phase hori- 
zon occurs. Because the observed radiation is thereby 
associated with this phase horizon, they claim it must be 
the Hawking effect. 

In the following, we shall examine this claim in more 
detail. More precisely, the question we are going to ex- 
amine is whether or not the existence of such phase hori- 
zons is enough for concluding the occurrence of Hawking 
radiation, and whether the existence of such an horizon 
can be responsible for the radiation they see. Note that 
the system studied in [l7| is very different from most of 
the previously studied black hole analogues, since there 
is no low-frequency regime in which the phase and group 
velocities are equal and where the black hole horizon is 
unambiguously defined. A crucial feature of their sys- 
tem is that the dispersion relation of the electromagnetic 
field in the frequency regime of interest is far from trivial, 
and the group and phase velocities are different. Further- 
more, there is no group velocity horizon (in the frequency 
range of interest), i.e., the group velocity of the photons 
is always smaller than the pulse speed. 

Since it is hard to analyze a model in which the true 
dispersion relation of the silica glass [2^ is taken into 
account [2l|, we shall assume that what is crucial in the 
argumentation in p7| is the existence of a phase - but 
not group - velocity horizon, and look at a model which 
shares this feature of their system although both, the 
cause of the dispersion relation, and the exact analytic 
form will be different. Because their claim (that they 
have seen Hawking radiation) rests crucially only on the 
existence of a phase horizon, the fact that this phase 
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horizon is due to a different mechanism in our model from 
theirs should not make any difference if their argument 
was correct. 

We shall use two models which behave very similarly to 
their model in the frequency/ wave-number regime of in- 
terest. These models will have phase horizons in exactly 
their sense, with the same characteristic parameters ~ but 
no group velocity horizons in that regime, just as theirs 
does not. The models clearly have an altered effective 
space-time geometry, and one can ask what the quantum 
radiation is from that space-time and whether it could 
be classed as "Hawking radiation" . 



II. EFFECTIVE METRIC 

Let us start with the macroscopic Maxwell equations 
in a medium whose index of refraction n is purely de- 
termined by the relative dielectric permittivity e — v?, 
i.e., the relative magnetic permeability is unity /i = 1. In 
natural units Eq = /io = cq = ft = 1, they read 

V -D = V • [n^E) =0, (1) 

V-B = 0, (2) 

Vxi? = V X B = dtD ^ dt{n^E) , (3) 

V X E ^ ~dtB . (4) 

Using the temporal gauge Aq = (which is not equivalent 
to the Coulomb gauge V • A = in this case), the two 
equations 1^ and can be satisfied automatically by 
introducing the usual vector potential via 



E = dtA, 
B = -V X A. 



(5) 
(6) 



The remaining two Maxwell equations, i.e.. Coulomb's 
law ([T]) and Ampere's law ([3]), then become 

V • (n^E) = V • (n'^dtA) = , (7) 
dt{n^E) = dtiri^dtA) = V x B = - V x (V x A) 



V'^A-V{V-A). 



(8) 



Let us now chose a specific polarization. For simplicity, 
we assume that n is a function of t and z only 



n = n{t, z) . 



(9) 



Thus the kx and ky, the eigenvalues of the operators idx 
and idy, are constants - which allows us to use the separa- 
tion ansatz A{t, r ) = A(t, z) exp{ikxX+ikyy} . Given this 
ansatz, we can always rotate the system so that ky = 0, 
and the solution A{t, f) is independent of y. We now 
choose the polarization with A^ = A^ ^ 0. Since A{t, f) 
does not depend on y, this automatically satisfies the 
Coulomb law ([7]) and Ampere's law ([8]) becomes 



where A = Acy and V has only x and z derivatives. 

For the other polarization, it is more convenient to 
employ the dual potential A such that 



D = n^E = V X A, 
H = B = dtA. 



(11) 
(12) 



This automatically obeys Coulomb's law ([T]) and 
Ampere's law ([3]). The remaining Maxwell equations read 



V • B = V • (<9f A) = , 



1 



(13) 



V X E = Vx|^— VxAj= -dtB = --d^A . (14) 

We see that, for a purely dielectric medium without mag- 
netic response, the description in terms of the dual po- 
tential A is actually simpler since it obeys the analogue 
of the Coulomb gauge V • A = in view of Eq. (fT3)) . In 
complete analogy to the previous case, we select a so- 
lution A which is independent of y and, consistent with 
V • A = 0, we take A^, = A^ = 0. With these simplifica- 
tions, Faraday's law ([H]) can be cast into the form 



(15) 



with A{t,r) = A{t, X, z)ey. We thus see that the two 
polarizations of the electromagnetic field do not obey the 
same wave equations pUj) and ([TSl) if n'^ actually depends 
on t or z (or both). 

The wave equations (fTO|) and p5|l can be written as the 
equations of motion of a scalar field in an effective 2-|-l 
dimensional (i, x, z) space-time. The effective metrics for 
the two polarizations are, cf. 



dsA 



1 



n^idx^ +dz^), 



(16) 
(17) 



dtin^dtA)=V'A, 



(10) 



where the first ds\ corresponds to the equation pO)) for 
A, and the second ds\ to the equation ([T5|) for A. The 
second effective metric is conformally related to the first 
so the light cones of both arc the same. Unlike the four- 
dimensional equations for the electromagnetic field, how- 
ever, these effective three-dimensional scalar equations 
are not conformally invariant and the two polarizations 
obey different equations of motion (for non-constant n). 



III. DISPERSION 

So far, we have described a dielectric medium without 
dispersion (where n depends on t and z, but not on w, 
for example). In a real medium, however, the electric dis- 
placement D(t,f) is not just given by the electric field 
E{t^ f) at that space-time point, but also depends on the 
electric field at earlier times. In a stationary medium. 



this can be represented via D{uj^ r) = e{uj, f)E{uj, f) af- 
ter a temporal Fourier transformation. For a dynamical 
medium, however, this is no longer possible in general be- 
cause of the time dependence of the dielectric constant. 

In order to avoid these difficulties, we shall mimic dis- 
persion by adding a mass term to the equations (|10p and 
P5|) of motion. As it is well known, a massive scalar 
field has a non-trivial dispersion relation, i.e., group and 
phase velocity depend on w. This simple model has the 
advantage that the generalization to curved space-times 
is straight forward. Equation ([TU]) is replaced by 



du 



(18) 



with the effective metric g'^'^ given by ([T5)) and analo- 
gously equation ([T5|) for A is obtained with the effective 
line element ([TT]). 

Furthermore, as in the experiment [l7| . we shall as- 
sume that the small perturbation of the dispersion re- 
lation (created via the non-linear Kerr effect) moves 
through the medium with a constant velocity v (deter- 
mined by the passage of an intense pulse of light through 
the medium). The two models we shall consider will differ 
in how that change is expressed. In the first model, the 
space-time dependence of the dispersion relation occurs 
through a variation of the refractive index 



n{z — vt) , 



(19) 



while the mass of the field is held constant. In the second 
model, the space-time dependence will come through a 
variation in the mass of the field 



m{z — vt) , 



(20) 



with the refractive index held constant. In the first 
model, the equation of motion (ITSl) has the form 

[dtn'^iz - vt)dt -dl- dl + ^^(z - vt)mf\ A = Q. (21) 

In the second model, the refractive index is assumed to be 
constant n = uq, but the mass m changes as a function 
of space-time 
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1 



A = 0. 



(22) 



In order to make this model correspond most closely to 
the experiment [l3|, we choose 



lim n 



lim m 



no 



1.4595, 

0.208 
/im^ 



2tt 



13.46/im / 



(23) 
(24) 



These factors were chosen so as to make the asymptotic 
effective refractive index noff(a;), given by the inverse of 
the phase velocity of the waves 



as close as possible to the actual refractive index mea- 
sured in the experiment [l3l. In particular the above 
values of no and mo were chosen to fit the Silica glass 
refractive index at 700 and 1100 nm. Note that A is the 
wavelength of the light in vacuum, not in the medium, 
i.e., it is really shorthand for 27rco/w. 

The wave equations for the other polarization read 



1 



n?{z-vt) ^ n^{z-vt) 

for the first model, and, for the second one. 



A = 0, (26) 



d 



1 



{dt+dt) + 



,2n , m^iz-vt) 



A = 0. 



(27) 



Compared to the wave equations (PT|) and (P^. the dis- 
persion relation is the same with the mass term divided 
by n^. Thus, in order to fit the dispersion relation of 
the silica glass, we can use the same parameters for 
no = 1.459 and take = m\rL^ = 0.943/^m^. 

Of course, Eq. (jl8p is simply the equation of motion 
of a massive scalar field in a background effective metric. 
As is well known, the phase velocity of a massive field is 
everywhere greater than the velocity of light (velocity of 
the characteristics of the metric) , and the group velocity 
is everywhere lower, but the product is always equal to 
the velocity of light (in this case c = 1/no). 

We are not saying that such an effective mass is the 
actual mechanism which determines the refractive index 
in the medium. It is clearly not. Rather it is the reso- 
nances of various electronic transitions which determine 
the refractive index of the system. However, our model 
mimics those features of the medium which the authors 
of [TtI claim to be important to their explanation of the 
phenomenon they observe. 

In particular, outside the band of wavelengths from 
700 nm to 1100 nm this effective refractive index deviates 
significantly from that of Silica glass, and even within 
this band the second derivative of the refractive index 
with respect to the wavelength has the wrong sign, see 
Fig. [H However, no- where in the analysis of [I4I does 
the behavior of the refractive index outside this band, or 
the sign of the second derivative of the refractive index 
play a role. It they were important, that dependence 
would in itself case serious doubt on the attribution of 
the radiation measured in [l3] to the Hawking effect. 



IV. FIRST MODEL 

In the experiment [l^, the velocity of the pulse v 
is larger than the velocity c = 1/no by a factor of 
1.4590/1.4533 ~ 1.004. Thus, we may apply a coordinate 
transformation (an effective Lorcntz boost) after which 
the pulse is instantaneous. To this end, we define a new 
time coordinate {v is assumed to be constant) 
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FIG. 1: Wavelength dependence of the effective refractive in- 
dices Tiph = c/i)ph, Tigr — c/vgT, and no — c/ y^iiphVgr for sihca 
glass (black) and for the massive approximation (green) which 
tries to fit the phase index of silica glass in the regime of inter- 
est by giving the field an effective mass. While this matches 
the phase index not too badly, it is worse for the group and 
geometric mean. 



such that n depends on t only n = n(T). In order to 
obtain a diagonal metric, we also introduce a new spatial 
coordinate 



an effective anisotropic "Hubble parameter" given by 
1 d\n{gpp) n'^v^ c?ln(n) 



2 dr 
dln(n) 
dr ' 



1 dr 



(31) 
(32) 



Note that we arc measuring the "Hubble parameter" us- 
ing the time dr which is related to the laboratory time dt, 
rather than the proper time ds of an observer at "rest" 
{p ~ const). Since n^v^ is quite close to unity, the two 
"Hubble parameters" are very different ^ Ux- 
The equation of motion (for one polarization) reads 



dr 



dr- 



2 2 

11 m 



A = (33) 



Since dp and dx correspond to Killing vectors, i.e., sym- 
metries of the metric, we may employ the separation 
ansatz A{t, p, x) = A{t) expjiKp -I- ik^x} which gives 



dr 



2 2 



1 



dr 



1 



2 I I 2 I 2 2 



A = (34) 



In order to study the solutions of this equation, it is useful 
to introduce yet another time coordinate via 



^2^2 _ I 



(35) 



In terms of this time coordinate, each mode {n^kx) cor- 
responds to a harmonic oscillator 



(- 



rf2 



f7^(T) M =0, 



with a time-dependent potential 



(36) 



p = z — I dr 



v^n'^{T) — 1 



(29) 



With the parameters from the experiment [T3| , the factor 
v'^n^ir) is always greater than unity - even for the max- 
imum value of 6n: for n = hq ~ 1.459 and v = c/1.4533 
as in [l3|, it is 1.008, whereas 5n « 0.001. 

Thus the above coordinate transformation is non- 
singular and results is a regular metric 



dsi 



2 2 



■dr' 



■ dp^ - n^dx^ , (30) 



as well as ds\ 



ds\/n^ 



j^l II for the other polarization. Re- 
membering that n depends on r only n = n{T), we 
see that this metric is purely time-dependent and corre- 
sponds to an expanding or contracting Universe - with- 
out any black-hole horizon. Hence, at best one could ar- 
gue that, if one saw particle creation in this system, one 
would have an analog to cosmological particle creation 
(see, e.g., [H, [l^]), not black hole evaporation. There is 



n^K^ + [kl + n^m^) 



(37) 



Again, since n^v^ is quite close to unity, the first term 
n^K^ will dominate unless k is very small. Before the 
pulse arrives, we have a constant potential 



if 



2 2 I /; 2 , 2 2 



1 



(38) 



and exp{— iilo^ + inp + ik^x} is a solution. Translation 
to laboratory coordinates exp{—iiot + ik^z + ik^x} yields 
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rig + nlv^K 



9 9 



(39) 
(40) 



After the pulse, the potential is constant ft^ = JIq 
again. However, a solution which initially behaves as 
A = expl—ifloT + inp + ikxx} will after the pulse be a 
mixture of positive and negative frequencies in general 
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with the Bogoliubov coefficients Q!„_fc^ and PK.,k^- The 
probabihty for particle creation (per mode k, k^) is given 
by 1/3^ Since the variation Jfi of the potential is small, 
we may use perturbation theory to estimate the expected 
number of created photons per mode via [1^, [2^ 



j dT5n{T) cxp{2inoT}. 



(42) 



From Eq. ([57]) . wc obtain 5fl ^ 5n. Thus, in view of 
the smallness oi Sn w 10""^, we can only hope for signif- 
icant particle creation if the T-integration compensates 
this small number. If we assume the time-dependence 
of 6rt{T) to be pulse-like (e.g., Gaussian), this is only 
possible if Qq is small enough, i.e., if Sfl constitutes a 
sizable fraction of flo- As discussed after Eq. ([ST]), this 
requires a small k. Thus, let us take k = in the follow- 
ing. Note that k = does not imply fc^ = according to 
Eq. (|40)) . Setting k = and = 0, we obtain the mini- 



mum value of Hq to n^m^ 
which yields 



1/u 



1 



"0 



1 



277 



1 1 [m\ 70 nm. 



©(lO^^Hz) 



(43) 



In the experiment [13], '^(t) changes on time-scales Ar 
of order lO^m. However, after the transformation ([55]). 
we get a rescaled period AT of order 2500//m. Conse- 
quently, since f2™™AT 1 is very large, the particle 
creation amplitude (|42|) is exponentially suppressed by 
many orders of magnitude. 

On the other hand, if n(j) was to change on optical 
time-scales Ar of order /im, i.e., femtosecond, we would 
get a rescaled period AT of a few hundred femtoseconds. 
Compared to the above value of Slg"", this would just be 
about the same order of magnitude - i.e., the particle cre- 
ation probability would not be too strongly suppressed. 
Nevertheless, if we estimate the order of magnitude of 



27r 



1300/im ' 



(44) 



we find that the probability j, | could be at the level of 
one percent at most. Still, one might hope to observe the 
analogue of cosmolo gica l particle creation with a set-up 
as in the experiment [l7|, provided that n(r) changes fast 
enough. However, if we insert a transversal wavenumber 
in the optical regime, such as 



k..^O\^-^ 



(45) 



we increase J7 by a factor of around six and thus dimin- 
ish the particle creation probability drastically (e.g., for 
a Gaussian pulse, several orders of magnitude). Conse- 
quently, the photons that are created would be emitted 
predominantly in forward direction and not to the side. 

The analysis for the other polarization is very simi- 
lar. Apart from the conformal factor the effective 



metric ds\ is the same as in Eq. (pO)) . Thus the wave 
equation reads 

n'^v^ - 1 

Or TT^r- Or T 



1 P V? ^ n4 



A = (46) 



Again, we can introduce a new time coordinate 



T= dT 



1 



(47) 



and thereby map each mode to a harmonic oscillator with 
a time-dependent potential. Since the change in n is 
quite small Sn « lO"'^ and mainly the relative change 
of n^v'^ — 1 matters, the additional factors of do not 
qualitatively affect our main conclusions. 



V. SECOND MODEL 

For the above model we have assumed that the mass 
of the field is constant in space and time, and that the 
change in the effective index of refraction n(uj) = \k\/u) 
is caused by a change in the effective velocity of light 
c = co/n. We could also choose to have the change in 
the mass generate the change in the effective index of 
refraction. In this case that change would not be inde- 
pendent of A but would scale as A^. However we also 
know that the change in refractive index of silica glass 
due to the Kerr non-linearity will not be independent of 
A, and nothing in the analysis of [l7| was sensitive to any 
such dependence of the Kerr non-linearity on frequency. 
Thus our choosing such a variation should not take this 
model outside the range of conditions assumed in [l7j . 

In the laboratory frame we have 



n{uj) = — = , 



so a change Sm in the effective mass gives 



Sn 



km 



: 6m 



km 



■ 6m . 



(48) 



(49) 



Reproducing the variation (5n « 10 at optical frequen- 
cies uj = 0{2Tr/ fim) and wavenumbcrs k requires 
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100/im 



(50) 



for one polarization A, for the other A, it is about a 
factor of two larger. The remaining analysis is analogous 
to the previous section. Again, particle creation is most 
pronounced for k = k^ = 0. The variation 6m induces a 
small change of the effective potential 



6n 



27r 



560/^m ' 



(51) 



for one polarization A and 6n « 27r/(280/im) for the 
other A. Since the minimum value il™'" is basically the 
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same as in the previous section, we again find that the 
particle creation probabihty is negUgible for the param- 
eter At — O{10fim) used in the experiment [o]- A 
measurable photon emission requires changing TO(r) on 
optical time-scales (or faster). In this case, the particle 
creation probability would be a bit higher than in the 
first model, but still not more that a few percent. 

But clearly any particles created simply by a change 
6m in the effective mass of the the propagating scalar 
particles is not the Hawking effect. This is related to our 
previous point, i.e., the absence of an effective black hole 
horizon. Instead, the photon creation mechanism is more 
similar to cosmological particle production. This effect 
can often be mapped to a time-dependent effective mass. 

On the basis of the experiment there is nothing 
to choose between either of the above models. But both 
of those models do not produce any particles due to an 
analogue of the Hawking effect. Furthermore, the num- 
ber of particles created, and especially the number of 
particles created in a direction perpendicular (in the lab- 
oratory frame) to the propagation of the pulse (i.e., with 
kx 3> fc^ ) ar e far too small to account for the observa- 
tions in [131 ■ In addition, the created photons do not 
necessarily occur in the frequency band the authors of 
1171 observe. 



the modes could always be arranged (by suitable rota- 
tion), this is not true in this case. Any ip dependence can- 
not be eliminated by a coordinate transformation. How- 
ever, the dominant photon creation would typically occur 
for the if independent modes with only minor amounts 
from the modes with higher angular momentum. Thus, 
we focus on the ip independent and dominant solutions 
of Eq. (|54p in the following. 

We still have the freedom of selecting a polarization. 
If we choose = Q (analogous to a TM mode), we get 
Az = and from V • A 0, we find A,, = 0. We are 
finally left with A^(t, r, z) and after adding the effective 
mass term, Eq. ([54)) becomes 



d? - rdr^dr - sAs. + ^ I A, 



0, 



(55) 



where we have added the effective mass term in order to 
model dispersion. 

The other polarization with = (analogous to a TE 
mode) can be described in complete analogy in terms of 
the vector potential with = 0. Again focussing on 
the dominant ip independent modes, V • {n^df A) ~ 0, 
i.e.. Coulomb's law ([7|), implies = 0. We arc left with 
A^{t, r, z) and Ampere's law ([5]) becomes 



VI. FILAMENTS 

So far, we have assumed planar symmetry, i.e., n and 
m did only depend on t and z, but not on x or y. How- 
ever, in the experiment [13, the travelling pulse is not a 
plane fronted wave, but is a thin filament. I.e., the pulse 
of changing n and m docs not only depend ont — z/v but 
also on X and y. In general the reduction of the Maxwell 
equations to an effective scalar field equation is much 
more difficult than in the planar case above. To sim- 
plify the analysis, we assume rotational symmetry which 
should be a reasonably good approximation. Thus, in 
cylindrical coordinates with r = \J x^ + U'^, we have 



dtn^dt - rdr-dr - dt + n^m' 
r 



A^^O. 



(56) 



Note that we are using the co- and contra-variant nota- 
tion Af = A^jr^ and A*^ = A^^/r^ where A^ and A^ 
must go to zero as or faster in order that E and B be 
regular at r = 0. 

The above equations (|55|) and ([56]) can again be in- 
terpreted as scalar wave equations in 3+1 dimensional 
curved space-times with the effective metrics 



d^ 



n 



2 I 2 

n ar 
dr^_ 



2j 2 

n dz 



d^ 



dz^ 



(57) 
(58) 



n 
m 



n{t,r, z) ^ n{t — z/v,r) . (52) 
m{t,r,z) ~ m(t ~ z/v,r) . (53) 



As the next step, we re-write Faraday's law ([Ti]) as 



provided that only ip independent solutions arc consid- 
ered. The weird 1/r^ dependence of the the angular d(p'^ 
term stems from the re-interpretation of the double vec- 
tor product as in Eq. ([54| as a scalar Laplace operator. 
Again, making a suitable coordinate transformation 



X I —V X A 



-V-A. 



(54) 



where I? is a self-adjoint operator acting on the Hilbert 
space of all transversal V • A = vector- valued functions 
A. Since n does not depend on ip, this operator 2? com- 
mutes with the generator Lz = id^ of rotations around 
the z-axis. Thus, we may classify the solutions of Eq. ([M)) 
in terms of id^ eigenmodes. 

In contrast to the planar case, where one could argue 
due to the rotational symmetry, that y-indepcndencc of 



\Jnlv'^ - 1 



t - 



P 



TIqVZ 



1 



(59) 
(60) 



we may diagonalize the unperturbed part of the metric. 
The resulting wave equation depends on the way how we 
model the change of the medium induced by the pulse - 
as a variation of the refractive index (first model) or via a 
change of the effective mass term (second model) . Since 
the latter choice gave us a slightly larger probability for 
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particle creation in the previous Section and is technically 
simpler than the former, we use the second model in the 
following. From the transformed metric 



ds\ = dr^ - dp^ - nldr^ - , 



(61) 



we obtain the wave equation 



1 



dt-d'p-—dr-dr + m'iT,r)\A^^O. (62) 



In analogy to the previous Sections, we use perturbation 
theory m = mo + dm where the unperturbed modes can 
be obtained by the separation ansatz 

An,.MAr, P, r) = Afn..,k^e~'''^ e'^PrMKr) , (63) 

with the Besscl function Ji. The Normalization factor 
J^ii,K.kr is chosen according to the pscudo norm 

(A\A) = - / dpdr—(A*drA-AdrA*) , (64) 
2 J r 

discussed in the appendix, which gives 



The frequency fl is given by 



rig 



(65) 



(66) 



Within perturbation theory (see appendix), the first- 
order amplitude -4^;"^^'''^, for creating a photon with 
il, K, kr from an initial quantum vacuum fluctuation with 
17', k', k'r is given by the overlap integral 

-^o-'k^I-' [ drdrdp An',K',A;;^o,K,^v ■ (67) 



As the particle creation process is basically multi-mode 
squeezing, the same ■^n'^K'^fe' yields the amplitude for 
creating a pair of photons with the quantum numbers 
fl, K, kr and ft', k', k'^, respectively. 

Since 6m does not depend on p, the p integration yields 



Thus, the particles are emitted in opposite p 



directions. Of course, this is strictly true only for an 
eternally propagating pulse Sm = 5m(t — z/v,r), a finite 
life-time would also induce a weak dependence on p. 

In the experiment [l7l | , 5m varies on time scales of the 
order of 10/xm. Again, the rate of change becomes even 
slower At = O{2Q0pm) after the transformation ((59)) 
and ()60|) to the new coordinates t and p. In view of the 
value for toq = 27r/(13.46//m), the frequencies Q, + U,' oi 
the photons to be emitted are far to large in compari- 
son to At and hence the integral over r is exponentially 
suppressed by many orders of magnitude. 

Similar to the planar case discussed before, we can only 
hope for a measurable particle creation probability if 5m 



varies on optical time scales (or faster). However, even 
in this case, significant particle creation is only possible 
for small enough k, and kr . Consequently, the transversal 
wavenumber kr should be far below the optical regime - 
i.e., optical photons are again predominantly emitted in 
forward direction. Since the transversal extent A?' of the 
pulse is a few microns at most - and thus much smaller 
than 27r/fcr ~ we may approximate the Bcssel functions 
describing the radial dependence via Ji{krr) sa krr/2. 
After that, the remaining power law dependence on kr 
and k'r can be pulled out of the integral and the total 
amplitude ^qJ'^^^^j,/ factorizcs. The dependence on k and 
k' is given by the 5{k -\~ k!) term mentioned above. The 
remaining term depending on -f- 17' reads 



^^g-»(n+a')r / drr^5m{T,r). 



(68) 



Thus, in comparison with the planar case discussed in the 
previous Sections, we find that the pair creation probabil- 
ity (per unit time) is additionally suppressed by a factor 
of (Ar/Ar)4. 

The other polarization (A) behaves in a very similar 
way - the equations are modified by factors of np in var- 
ious places. However, that does not change the order of 
magnitude estimated above. 



VII. PULSE SUB-STRUCTURE 

So far, we have considered a rigidly moving pulse, i.e., 
n = n(t — z/v,r) and m = m{t — z/v,r). However, 
the pulse in the refractive index created by the incoming 
1.06/im radiation in [iTj is not a simple (exponential) 
pulse, but contains sub-structure. The index of refraction 
change 5n goes as the square of the incoming electric 
field strength 5n oc , and since the incoming radiation 
is linearly, not circularly, polarized, the electric field E 
oscillates at the frequency of the incoming radiation, and 
the square then oscillates at twice that frequency. 
Taking this into account, a more realistic model of 5n 
will then have the form (and similarly for 5m) 

5n{t,z,r) = (5no(r) cxp | - , ^^^} > x 



X cos Win < t 



2(Ai)2 

z 

Vph 



(69) 



Here is the frequency of the incoming radiation 
(which generates the pulse) and its phase velocity 
- corrected by a factor of 1/ cos 8, where w 6.5° is the 
cone angle of the Besscl pulse. 

The oscillating pulse sub-structure given in the second 
line has two important consequences. First, in contrast 
to the pulse envelope, which is slowly varying on a time 
scale At > O(lO^m), the oscillating term is changing on 
optical time-scales (with 2a;i„), i.e., much faster. Second, 
since the pulse speed v « cq/ 1.453 and the phase velocity 
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■^ph ~ Co/1.44 do not coincide, the perturbation does not 
just depend out (x {t— z/v), but also on p. Even though 
the difference between v and iiph is only on the percent 
level, this deviation is enhanced by the coordinate trans- 
formation ((59)) and (|60p to the r, p coordinates 



Sn(T,p,r) — 6no{r)c:xjp 



2(novAt)- 



> X 



X COS^ (Wi-T + UJpp) , 



(70) 



with the frequencies 



1 



Win 



now 

- V/Vph 



1 - w/wph 



27r 

5/im ' 



1 



27r 
10/im 



(71) 



VIII. CONCLUSIONS 

There are a number of conclusions one can draw from 
our analysis for the interpretation of the experiment jl7| . 

First, whatever is causing the particle production due 
to the intense pulse, it is certainly not Hawking radia- 
tion caused by the "phase horizon" . The effective metric 
does not correspond to a black-hole horizon but has more 
similarities to a cosmological set ting - similar to a purely 
super- luminal pulse, see also [i^,!!^- Simulating a black- 
hole horizon would require pulses with other parameters, 
for example a stronger non-linearity Sn^ such that the 
effective metric in (pO| is no longer regular. 

Second, even the interpretation as cosmological parti- 
cle creation yields negligible photon pair creation proba- 
bilities in the optical regime if we consider - as done in 
p7j - a rigidly moving pulse 6n{t — z/v,r) or 6n{t — z/v). 
After the effective Lorentz transformation, the rate of 
change of the pulse is just too slow. However, taking 
into account the oscillating sub-structure of the pulse 
(not just its envelope, as done in [l^l), one could get 
a non-negligible probability for particle creation. 

Third, in all of these scenarios, the creation of photons 
with transversal wave-numbers in the optical regime is 
strongly suppressed - almost all the photons are emitted 
roughly in longitudinal direction. Thus it is very hard to 
explain the photons emitted in perpendicular direction, 
as observed in [l3|. If one were to include the effects of 
defects in the dielectric medium, this could change. On 
the one hand, they could effectively scatter the co-moving 
photons out of the co-moving direction. However, since 
it is hard to imagine more than a tiny fraction of the 
photons being so scattered, the photon creation rate by 
the pulse would have to be even larger to create enough 
photons. On the other hand, these defects could also 
induce short-scale deviations from a rigidly moving pulse 
6n(t — z/v, r) with comparably high frequencies - similar 
to the pulse sub-structure discussed above - and thereby 
increase the particle creation rate. 



Fourth, a key observation is that the particle produc- 
tion rate is not simply determined by the structure of the 
phase horizon. It depends in detail on how that phase 
horizon is created. As shown in the case of the "slab" 
geometry for the pulse, the two models for the varia- 
tion of the material properties (first model versus second 
model) yield slightly different numbers of particles. This 
is another indication that this phase horizon is not an ap- 
propriate analog to the black hole horizons. In the latter 
case, the radiation (i.e., black hole evaporation) is solely 
a function of the structure of the horizon itself (i.e., the 
surface gravity). Here it depends in detail on how the 
horizon is modeled. 

Fifth, it is also of importance that the equations of 
motion of the two polarizations {A^ and A;^) are not the 
same. While the effective metrics in the two cases are 
conformally related to each other, the wave equations for 
the polarizations are not conformally invariant. Were the 
effective metrics real 3+1 dimensional metrics, and the 
fields were the electromagnetic fields in 34-1 dimensions, 
the equations would be conformally invariant. However 
here the effective metrics are in reduced dimensions, and 
the equations for the potentials are effective scalar field 
equations. Both destroy the conformal invariance of the 
equations of motion. Thus, the propagation and particle 
production rates for the two polarizations will differ. 

Finally, we find that more investigations are needed 
to determine the cause of the radiation observed in [l3| 
and its potential relation to quantum radiation phenom- 
ena such as cosmological particle creation. We cannot 
claim to have explained their observations. This analog 
model clearly needs further work. Our modelling of the 
dispersion relation of light in such a dielectric medium is 
at best only a crude model of what takes place. Mod- 
els which reproduce the dielectric behavior as encapsu- 
lated in the Sellmeier coefficients [1^ are possible, and 
are being studied. However, we do not believe that they 
will alter our conclusions. Even were they to increase to 
the particle emission rate in a perpendicular direction, 
they would simply emphasize the dependence of the phe- 
nomenon on the details of the model, rather than relying 
solely on the structure of the horizon. 
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Appendix 

We wish to find the lowest order amplitude for the 
particle creation rate caused by a small fluctuation of 
the mass Smir, r) in the equation 



dl-dl-^dr-dr+ml + 2modm ) * = . (72) 



In terms of the momentum 11 = 9,- 4*, the inner product 



Initially r — > -co, i.e., before the pulse arrives, the solu- 
tion (^P, n) is supposed to coincide with the unperturbed 
mode / = 0. Thus we have ao{T — > —oo) ~ 1 while 
all other a/ as well as all /?/ vanish for r — ^ — oo. In 
the final regime r — > oo, however, some of these modes 
will be excited (3i = 0{Sm) and a/^o = 0{Sm) due 
to the perturbation S'M. induced by the pulse. Inserting 
the expansion (j76p into the equations of motion (|75p and 
neglecting terms of 0{[Sm]'^), we find 



(«'i|*2) = ^y"drdp^(*jn2-nt*2) , (73) 

is conserved 9t(^'i|^2) = for any two solutions 4'i and 
^2- Introducing the matrices 



Mo 







1 



I?2 _ rng 



, SM = 





2'mo6m 



, (74) 



where ^ d'^ + ''r drr ^ dr is the spatial differential 



3p -r -i-o 

operator in Eq. ((72|). the equations of motion can be cast 
into the from 



E 



n/ 



drai 



n* 



drf3l 







(77) 



Projecting the the above equation onto the mode J via 
the inner product ([73]). we get 



d^l3j = -imo J drdp'^Sm^j^o + 0{[6mf) . (78) 
Integration over r yields Eq. (p7)) 



dr 



(Mo + 5M) 



(75) 



Now, let (^7,11/) be a complete set of solutions of the 
unperturbed wave equation (j72p with 6m = which have 
positive pseudo norm and are orthogonal (^'/l^'j) = Sij 
with respect to the inner product ([75]). and such that 
(\Ii^|v]/*) = 0. For simplicity, we choose the {^i,Ui) as 
eigenvectors of the matrix operator Mo with eigenvalues 
—irij. Since Mo is real, the complex conjugated solutions 
(^|,nj) then form a complete set of negative pseudo 
norm (^||^}) = ^Sjj. Thus we can expand the solution 
(^P,!!) of the full wave equation ([72l) with Srri ^ into 
these sets 



2n 

I3j(t -> oo) w -imo / dr dr dp — 5m ^-q . (79) 



This gives the amplitude for converting an initial quan- 
tum vacuum fluctuation in the mode 5*0 into a real par- 
ticle in the final mode '^j. The total probability for 
creating a particle in the final mode J is then given by 
the sum over all initial modes / 



I 



27r 

dr dr dp — (5m ^' / 
r 



(80) 



Pi{r) 



(76) For our system, this becomes 



PiJ = A^v1':k = --^sJ^ j drdrdpr5m{T,r)e-^^^e}^'>J,{Kr)e-^^'^e^-'''JMr) 

= + J dTe-(^^+"')-| drrSmiT,r)Mkrr)J^iKr). (81) 



Calculating the total probability Vn.K,kr then gives the 
usual (5(0) type singularity from the k' integral over [(5(k-|- 
k')]^ because of the assumed unboundedness of the pulse 
in the p direction. Removing it gives the particle emission 
rate per unit length of the pulse in the p direction. Note 



that 17 and Vl' are both positive, as we have chosen modes 
going as e"*^"^ as our set of positive pseudo norm modes. 
The Bessel functions Ji{krr) go as k^r /2 when small, and 
become oscillatory for large fc^r. 
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